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Abstract 

Cross sections are calculated for single and double photon exchange in ultraperipheral Pb+Pb 
collisions at the LHC. The particle production is simulated with the DPMJET event generator. 
Large cross sections are found for particle production around mid-rapidity making these processes 
an important background to hadronic nuclear interactions at both the trigger and analysis levels. 
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The strong electromagnetic fields present in high-energy nuclear collisions with impact 
parameters larger than the sum of the nuclear radii lead to large cross sections for a vari- 
ety of photonuclear processes. In these ultraperipheral collisions, there are no geometrical 
overlaps between the colliding nuclei, so purely hadronic interactions are suppressed. Parti- 
cle production in ultraperipheral collisions has been studied in Au+Au collisions at RHIC 
and also in heavy-ion interactions at lower energies. See Ref. [lj for a review. These stud- 
ies have so far been focused on the exclusive production of a single (vector) meson, e.g. 
Au + Au — > Au + Au + V with V = p° or J/^, or on two-photon production of dilepton 
pairs. The nuclei have remained intact or have only been excited to low energies by the 
exchange of an additional soft photon. 

In this paper, we consider particle production in a general photonuclear interaction, 
7 + A — > X in ultraperipheral collisions between two lead nuclei at the CERN Large Hadron 
Collider (LHC). We do the calculations for the collision energy ^/snn = 2.76 TeV which will 
be available during the first two heavy-ion runs in 2010 and 2011/12. As will be discussed 
further below, the relevant photon energies for particle production around mid-rapidity are 
between ~ 10 GeV and ~ 100 TeV in the rest-frame of the target nucleus. We consider 
the cases when one or two photons are exchanged. In the former case, the photon-emitting 
nucleus remains intact, A + A — > A + X, while in the latter case both nuclei disintegrate. 

The particle production in the photonuclear interactions is modeled using the DPMJET 
Monte Carlo event generator version 3.0 {2, 3|. 

We begin by discussing the photon flux associated with relativistic heavy ions and derive 
the relevant spectrum for single and double photon exchange. We then discuss particle 
production in these collisions and how the produced particles are distributed in phase space. 
We conclude by briefly discussing the experimental consequences of these processes in a 
typical high-energy collider experiment. 

The electromagnetic field of a relativistic particle can be treated as an equivalent flux 
of photons. This is the Weizsacker- Williams method. For collisions between relativistic 
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nuclei, the photon spectrum should be calculated in impact parameter space 
way, interactions with nuclear overlap, where the strong interaction dominates, can easily 
be excluded. The difference between this approach and using a nuclear form factor for 
calculating the photon spectrum will be discussed below. 

The photon spectrum, dn/dk, of one ion in an ultraperipheral collision is given by an 



integral over the impact parameter, b, 

The cut-off at small impact parameters ~ 2R is implemented using a smooth function for the 
probability of not having any hadronic interaction (1 — Phad(b))', Phadip) is calculated from 
a Glauber model. Since the photonuclear interaction probabilities are highest in grazing 
collisions with impact parameters b ~ 2R, it is important to take the hadronic interaction 
probability into account properly and not rely on a simple cutoff at b = 2R (see Fig. [T] 
below). The doubly differential photon spectrum d 3 n/dkdb 2 is 

d 3 n Z 2 a 1 2 2 . . 

dkd¥ = ^W x l[x} ( ) 

where x = bk/ r ) and 7 >> 1 is the Lorentz factor. 

The cross section for a photonuclear interaction with photons from a single beam is then 
given by 

roo 

(Ja+a^a+x = / —a lA (k)dk (3) 

The integral is cut off by the rapid decrease of the photon spectrum for k > 7/ R ~ 120 TeV 
in the rest frame of the target nucleus. 

By inserting the expression for dn/dk from Eq. (CQ) into Eq. ([3]) and changing the order 
of integration, the first order photonuclear interaction probability as a function of impact 
parameter can be obtained through differentiation: 

^ ^ db 2 Ik mi „ dkdb 2 ^ ^ ^ 

The photonuclear interaction probability as function of impact parameter is shown in 
Fig. [TJ for three different values of k m i n = 6 GeV, 1000 GeV, and 10 TeV. The minimum 
value, 6 GeV, is chosen as the lowest photon energy that can be handled by DPMJET. 
The effect of this minimum on the particle production in different rapidity intervals will 
be discussed further below. The photonuclear cross section, a 7J y(k), is the cross section for 
particle production calculated by DPMJET; it depends weakly on the photon energy, k, 
and increases monotonically from 8.1 to 12.9 mb over the relevant energy range (6 GeV - 
100 TeV). 

The total photon-proton cross section has not been measured above center-of-mass ener- 
gies of 200 GeV, and there is thus some uncertainty in the energy dependence at the highest 



energies The energy dependence of the photonuclear cross section should, however, be 
weaker because of shadowing, which implies that only the edges of the nuclei are affected 
by the increase. The contribution to the total cross section from photons with these high 
energies (k > 21 TeV) is furthermore rather low, so this uncertainty will not significantly 
affect the results presented here. 

The interaction probability has a maximum of ~15% for k m i n = 6 GeV in grazing colli- 
sions with impact parameter ^16 fm. For k m i n = 1000 GeV and 10 TeV, the corresponding 
maximum probabilities are 6% and 2%, respectively. These high probabilities make exchange 
of multiple photons in the same event likely. For cases where the collision time b/j is much 
smaller than the excitation time, 1/k, the sudden approximation applies, and the probabil- 



ity for multiple photon exchange factorizes 



6|. This method has been used previously to 



calculate the total cross sections for correlated forward-backward Coulomb dissociation Q. 
The cross section for exchanging one photon from each nucleus is thus given by 

a AA = 2tt / [Pi(fe)] 2 (l - P had (b))bdb (5) 

The photon energy spectrum for double excitation is obtained by differentiating Eq. (jSJ) 

with respect to the photon energies k\ and &2- The result can be written 

dcr „ f°° d 3 n .. . d 3 n n .,..,„ 

dk\dk2 Jo dk\db A dk2db z 

It is worth noting that there is a positive correlation between the two photon energies. The 

photon spectrum from a single nucleus under the requirement that the other nucleus emits 

a photon is obtained by integrating over either fci or k-i- The result is 

lH*f ^ ( ™ (1 - p - (6))M6 (7) 

The spectrum is thus weighted toward a smaller impact parameter by the photonuclear 
interaction probability P\ib\ resulting in a harder spectrum. This is in agreement with 
what was found for photonuclear vector meson production in coincidence with Coulomb 
break up [8|. 

Although the interaction probabilities for these photon energies are fairly high, they 
are sufficiently smaller than unity for the first order probabilities to be used. Unitarity 
restoration jfj has been estimated to reduce the cross sections with at most a few percent. 

A high energy photon may interact with a nucleon or nucleus in two different ways: it 
can appear as a bare photon, which interacts with a parton in the target, or it can first 



fluctuate to a qq pair, which then interacts with the target via the strong interaction, see 
e.g. Refs. 0, [hJ]. The importance of the hadronic component of the photon at high-energy 
electron-proton colliders has been pointed out earlier 9|. Since the quantum numbers of the 
photon are 1 , it tends to fluctuate to a virtual vector meson (vector meson dominance). 
The bulk of the photonuclear cross section and the soft particle production can be understood 
from resolved vector meson interactions. 

The DPMJET Monte Carlo event generator is based on the two-component dual 

parton model. The resolved part of the photon-hadron interaction is handled by the dual 
parton model through its implementation in PHOJET. The scaling from photon-hadron to 
photon-nucleus interactions, including shadowing, is treated within the framework of the 
Gribov-Glauber approximation. The high-p^ partonic processes are calculated from lowest 
order perturbative QCD. The model has been tested against data from fixed target 7r-nucleus, 
p+nucleus, and /i+nucleus experiments. It has been shown to reproduce the bulk particle 
production well [3]. It can handle low or intermediate photon virtualities, which is fine in 
this case, since the photons from the Weizsacker- Williams fields have very low virtualities 
Q 2 < (V-R) 2 ~ 10~ 3 GeV 2 . The fragmentation of the target nucleus is not considered in the 
current version of DPMJET, but knock-out protons from the target remnant are included. 

Samples of events have been generated with DPMJET with the photon spectrum de- 
scribed in the previous section and with different minimum photon energies. Six samples 
have been generated with 500,000 events per sample. Samples with single and double ex- 
citation have been generated with values of k min = 6, 1000, and 10,000 GeV. The cross 
sections for these samples are listed in Table [H The table also includes the corresponding 
cross sections when it is required that there should be at least one charged particle at central 
rapidities. 

The normalized charged particle pseudorapidity (77 = — ln(tan(6 l /2))) density distribu- 
tions for the samples are shown in Figure [2] (a) - (f). As expected, the particle production 
becomes more centered around mid-rapidity with increasing photon energy. The correspond- 
ing multiplicity distributions of charged particles with < 1 are shown in Figure [3] (a) - 

(f)- 

The cross sections in Table [I] are huge. The total photoproduction cross section for 
k > 6 GeV is about 3 times larger than the total hadronic cross section. If one requires 
at least one photoproduced charged particle within \r}\ < 1, the cross section is about 4 b 



or roughly 50% of the total hadronic cross section. As can be inferred from Fig. [3j the 
probability for having larger multiplicities around mid-rapidity is also very high. It is 
essential, especially for the higher values of k min , to calculate the photon spectrum in the 
impactparameter space. If a form factor (a convolution of a hard sphere and Yukawa form 
factor [8J) is used, the cross section for single excitation with k min = 6 GeV has been found 
to be 30% higher than in Table I. For k min = 1000 and 10,000 GeV, the differences are 73% 
and 185%, respectively. Using a nuclear form factor does not remove the contribution from 
events where the nuclei overlap. 

The transverse momenta of the photoproduced particles are typical for soft hadron- 
nucleus interactions. The pt distribution at mid-rapidity for single production with k m i n = 6 
GeV is shown in Figure H] (a). The distribution is essentially the same for all samples studied 
with a mean transverse momentum of < px >~ 450 MeV/c. 

The large cross sections and the characteristics of these interactions should make them 
an important background to hadronic nuclear interactions at colliders. Single excitations, 
which have the largest cross sections, are characterized by a strong asymmetry around 
mid-rapidity event-by-event. They can thus to some extent be rejected by requiring the 
presence of particles on either side of mid-rapidity. The cross sections for double excitations 
are lower, but these events have particles produced over a wider range of rapidities. They 
may therefore be harder to separate from hadronic interactions with low or intermediate 
multiplicities. 

The photon spectrum of course extends to, in principle, arbitrarily low values of k. The 
lowest k min = 6 GeV used here is therefore in that sense arbitrary. However, for particle pro- 
duction around mid-rapidity, photons with low energy do not contribute. This is illustrated 
in Figure H](b) and (c). Figure H](b) shows the ratio of the number of events with at least one 
charged particle within \r)\ <1 to all events as a function of photon energy. The distribution 
goes to zero around k ~ 150 GeV, and photons with energy lower than this do thus not 
contribute to the particle production within the two most central units of pseudorapidity. 
Figure S] (c) shows the same distribution for events with at least one charged particle within 
1 7/| <4.5. The ratio goes to zero above the lowest k m i n = 6 GeV used here. The photon 
energy range considered should thus give a complete description of photoproduction within 
the nine most central units of pseudorapidity, \rj\ <4.5. 

Unlike exclusive production of single mesons or dilepton pairs, inclusive photoproduction 
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in ultraperipheral collisions has so far attracted rather limited interest. Some early studies 
were done for RHIC Jjj], mainly to find the background rates for exclusive production. The 
cross sections are about an order of magnitude lower at RHIC energies (^/snn = 200 GeV) 
than at the LHC Inclusive photoproduction of mesons in high energy heavy-ion interactions 



were also studied in Ref. 



131 ] . The results are not directly comparable to the present work, 



because of a lower minimum photon energy, k m i n , and because only multiple photons hitting 
the same nucleus were considered. 

We have calculated the photon energy spectrum for single and double photonuclear ex- 
citations in Pb+Pb collisions at the LHC. Events have been generated with the DPMJET 
Monte Carlo according to these spectra. Large cross sections are found, also for particle 
production around mid-rapidity. 

The large cross sections and the non-zero probability for having charged particles pro- 
duced with intermediate or high around mid-rapidity make these events an important 
background to peripheral and semi-central Pb+Pb collisions at the LHC. They have to be 
taken into account in order to correctly determine what fraction of the total hadronic cross 
section a given event selection corresponds to. 

If the photonuclear events can be clearly separated from the hadronic events, e.g. using 
the rapidity gap between the photon-emitting nucleus and the produced particles, much 
interesting physics could be extracted from them. The cross section to produce a cc pair 
through 7+gluon fusion is, for example, nearly 1 b at the full LHC energy 14j. If these 
events can be measured, they would provide valuable information on the nuclear parton 
distribution functions. 

We thank Stefan Roesler (CERN) for providing the code for the DPMJET Monte Carlo 
and for help in getting it running. One of us (JN) would like to acknowledge useful discus- 
sions with Spencer Klein (Berkeley) and Anthony Baltz (Brookhaven) on multiple photon 
exchange in a single event. 
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TABLE I: Cross sections for particle production for single and double photonuclear excitations 
with different minimum photon energies, k m i n . The total cross sections and the cross sections 
for having at least one charged particle within \rj\ < 1 are shown. The cross section for single 
excitations take into account that both nuclei can act as photon-emitter and target. 
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FIG. 1: Probability for having a photonuclear interaction and no hadronic interaction as a function 
of impact parameter for three different minimum photon energies, k m i n . 
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FIG. 2: Pseudorapidity distributions of charged particles for different minimum photon energies 
for single (left) and double (right) excitations with different cut-off energies, k m i n . For single 
excitations, the solid histograms are for photoproduction off the nucleus with negative rapidity 
and the dashed histogram corresponds to production off the nucleus with positive rapidity. The 
peaks around rj = ± 9 are from knock-out protons from the target nucleus. 
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FIG. 3: Multiplicity distributions of charged particles at mid-rapity \rj\ < 1 for single (left) and 
double (right) excitations with different cut-off energies, k m i n . 
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FIG. 4: (a) Transverse momentum distribution of charged particles within \r]\ < 1. (b) Fraction of 
events with at least one charged particle within \r/\ < 1 as function of photon energy, (c) Same as 
(b) but for \r}\ < 4.5. 
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